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ABSTRACT 

The frequency of heating events in the corona is an important constraint on the coronal heating mechanisms. 
Observations indicate that the intensities and velocities measured in active region cores are effectively steady, suggesting 
that heating events occur rapidly enough to keep high temperature active region loops close to equilibrium. In this 
paper, we couple observations of Active Region 10955 made with XRT and EIS on Hinode to test a simple steady 
heating model. First we calculate the differential emission measure of the apex region of the loops in the active region 
core. We find the DEM to be broad and peaked around 3MK. We then determine the densities in the corresponding 
footpoint regions. Using potential field extrapolations to approximate the loop lengths and the density-sensitive line 
ratios to infer the magnitude of the heating, we build a steady heating model for the active region core and find that we 
can match the general properties of the observed DEM for the temperature range of 6.3 < Log T < 6.7. This model, 
for the first time, accounts for the base pressure, loop length, and distribution of apex temperatures of the core loops. 
We find that the density-sensitive spectral line intensities and the bulk of the hot emission in the active region core 
are consistent with steady heating. We also find, however, that the steady heating model cannot address the emission 
observed at lower temperatures. This emission may be due to foreground or background structures, or may indicate 
that the heating in the core is more complicated. Different heating scenarios must be tested to determine if they have 
the same level of agreement. 
Subject headings: Sun: corona 



1. Introduction 

A major problem in coronal physics is to determine the mech- 
anisms that transfer and dissipate energy that heat the million- 
degree corona. Regardless of the transfer or dissipation mech- 
anism, the predicted timescale for energy release for all mech- 
anism s is finite, i.e., a single heating event is relatively short- 
lived (|Klimchukll2006l l . The frequency of heating events occur- 
ring on a single strand in the corona, however, could distinguish 
between the different theoretical mechanisms. (Here, we use 
the term "strand" to refer to the fundamental flux tube in the 



intensities and velocities. If the heating events are infrequent 
(i.e., the time between heating events is much longer than the 
cooling time of the plasma), then the resulting temperature 
and density of the plasma along the strand will evolve; such 
a heating scenario is oft en termed "nanoflare heating" (e.g., 
ICardU Klimchuk|[l997l ). though it could be representative of 
different types of mechanisms, not simply Parker's canonical 
nanoflare theory (Parker 1972). The evolution of the plasma 
may or may not be apparent in loop observations depending on 
the number of strands the make up the loop and the relative 
timing of their heating events. In the case of a loop formed of 



Winebarger et al. 

structures (iReald I2010I I. So-called "warm," 1 MK loops 
that form a bright arcade in EUV images are evolving and 
their properties are consistent with l ow-frequency heating in 
the form of short- nanofl are storms (Winebarger et al. 
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Warren et al. 2003; Winebarger V Varren 2005; Ugarte-Urra et al 
20061 . l2009l : lMulu-Moore et al.ll201lh . The intensities of the hot 



(> 2 MK) loops that make up the active region cores, however, 
appear to be steady over many hours of observation (e.g., 
I Warren et al. 20 IG. 120 111 ). This steadiness is apparent in both 
the X-ray intensity over the neutral line and in the footpoints of 
the h ot loops observed in th e EUV and commonly called moss 
fe.g.. ^Antiochos et al.l l2003l V Additionally, the velocities and 
non-thermal velociti es in the moss are also steady for hours of 
observations ( Brook s^ WarrenI 12009V The steadiness of the 
intensities and velocities has prompted many to hypothesize 
the heating in the core regions is effectively steady. However, 
different heating scenarios, such as the long nanoflare storm 
scenario described above, also predict steady intensities and 
velocities. Additional comparisons between observations and 
simulations are necessary to discern whether the active region 
core is consistent with effectively steady heating. I Warren et al.l 
(HoTT) demonstrated that the properties of the plasma along 
strands suffering effectively steady (high-frequency) heating 
can be well approximated by the steady heating solutions to 
the hydrodynamic equations. In this paper, we will determine 
whether observations of the loops that form the active region 
core are consistent with the predictions of steady solutions to 
the hydrodynamic equations. 

If the density and temperature in a well-isolated and re- 
solved loop was known, it would be possible to test whether 
the loop was heated steadily or not using the so-called Rosner- 
Tucker-Vaiana (R TV) scaling laws (|Rosner et al also see 

ISerio et aP Il98lh . These scaling laws define a relationship 
between the apex temperature, Tapex, base pressure, po, and 
loop half-length, L, for steady uniform heating, i.e., Tapex = 
1.4 X lO^(poT)^/^. However, individual loops are very difficult 
to isolate in active region cores, and loops that appear to be 
resolved may themselves be formed of sub-resolution strands, 
hence such a direct comparison with model predictions is not 
possible. 

Another method to test whether a heating scenario is viable 
is to forward model the active region or full Sun as ensembles 
of strands with either a steady or infrequent heating model and 
then compare simulated images from the model to the observed 
images . Several studies have been completed using this tech- 
nique ([Schri iver et al.l I2OO4I : IWarren Winebarged[2006l . l2007l : 
iLundquist e t al. 20081). In these studies, the field lines from the 
potential magnetic field were used to approximate the global 
or active region structure. The heating rate was assumed to 
be a function of the magnetic field strength and loop length. 
Several different parameters were varied, such as the functional 
form of the heating, the degree in which the loops expanded. 



emission was too bright. 

In an attempt to simultaneously match both the EUV 
moss emission and the apex X-ray emission in two filters, 
IWinebarger et aL 1 (|2008h modeled an active region core as an 
ensemble of strands. Instead of assuming a functional form of 
the heating, however, they used the TRACE 171 A intensity in 
the moss region to limit the steady heating rate. The TRACE 
171 A intensit y is pressure sensitive given an unknown filling 
factor ((Martens et ahlbOOOl V In this study, a filling factor was 
assumed for the entire moss region, then the pressure was de- 
rived from the TRACE 171 A image. For each strand, a heating 
rate was determined that would match the pressure; the ex- 
pected intensities in the X-ray filter data was then calculated. 
This process was completed for several different assumed filling 
factors. The best agreement for the X-ray filter intensities was 
found for loops that expand inversely with the magnetic field 
strength and for a filling factor of 8%. This study was the first 
successful steady heating model to account for both the foot- 
point EUV emission and apex X-ray emission. The weakness 
of this study was the forced assumption of a single filling factor 
for the entire moss region. 

In this paper, we will determine whether observations from 
Hinode XRT and EIS agree with a steady heating model of 
the active region core. We build upon the previous full models 
of active region cores with two key improvements. First, in- 
stead of comparing model intensities with observed intensities 
in one or two filter images, we compare a differential emis- 
sion measure (DEM) from the model with a DEM calculated 
from a combination of XRT filter and EIS spectral line intensi- 
ties. Second, we use two density sensitive spectral lines (Fexii 
186.880/195.119 A) in the moss region to simultaneously con- 
strain both the pressure in the core loops and the filling factor. 
As with the previous studies, we use potential field extrapola- 
tions to approximate the loop length and geometry, assume the 
loops expand with height and assume the energy is deposited 
uniformly along the loop. We determine that the DEM from 
the steady heating model well matches the DEM from obser- 
vations around the peak temperature in the active region core 
(6.3 < Log T < 6.7), as well as the spectral line intensities in 
both the Fexii 186.880 and 195.119 A lines in the footpoints 
of the active region. However, the steady heating model does 
not reproduce the observed DEM in the "warm" temperature 
range (6.1 < Log T < 6.3). This warm emission could orig- 
inate from the overlying arcade and not from the core loops, 
or it could indicate that the heating in the active region core 
is more complicated. These results provide, for the first time, 
evidence that the high temperature loops in the core of an ac- 
tive region and their footpoint pressure can be well represented 
by high-frequency heating. Different heating scenarios must be 
tested to determine if they have the same level of agreement. 

2. Data 
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ters that can be used alone or in combination. The XRT data 
used for this analysis were taken on 2007 May 13 from 16:17 - 
16:30 UT. The data set is comprised of full Sun, full resolution 
images in multiple filters and filter combinations. The data pro- 
cessing included dark-frame subtraction, vignetting correction, 
and high-frequency noise removal using the standard xrt_prep 
routine available from Solar Soft. Spacecraft jitter was removed 
using xrt_ jitter, and long and short exposures (see Table [U 
were co-aligned and combined to increase the image dynamic 
range. Additional Fourier filtering was done to the data from 
the thickest channel, Be_thick, to remove low-level, residual, 
longer-wavelength noise patterns. We have used updated filter 
calibrations (Narukage et aLllioill ) and accounted for a 1640- 
A thickness of the time-dependent contamination layer on the 
CCD. (We used the "best-model" contaminant composition of 
diethylhexyl phthalate, as provided by the XRT calibration 
software, but we note that the true composition is as of yet 
undetermined.) 

The EIS instrument on Hinode is a high spatial and spectral 
resolution imaging spectrograph. EIS observes two wavelength 
ranges, 171-212 A and 245-291 A, with a spectral resolution of 
about 22 mA and a spatial resolution of about l" per pixel. 
There are 1'' and 2'' slits as well as 40^^ and 266^^ slots avail- 
able. Solar images can be made using one of the slots or by 
stepping one of the slits over a region of the Sun. Telemetry 
constraints generally limit the spatial and spectral coverage of 
an observation. See Culhane et al.l (|2007l ) and lKorendvke et al.l 
dioos) for more details on the EIS instrument. 

In this analysis we consider EIS observations taken from 
14:25:57 to 15:15:10 UT on 2007 May 13. For these obser- 
vations the 1" slit was stepped across the active region and a 
10 s exposure was taken at each position. The area observed 
was 256^^ X 256^^ A total of 16 data windows, some containing 
multiple emission lines, were downloaded from the spacecraft. 
The raw data was processed to remove the CCD pedestal, dark 
current, and spurious intensities from warm pixels. The pre- 
flight calibration was also applied to the data. Before the line 
profiles are fit it is necessary to correct for an oscillation in the 
EIS wavelength dispersion. We do this by assuming that there 
are no net velocity shifts along the bottom 30 pixels of the 
slit fsee lBrown et al1l2007l for details). For each emission line 
of interest the best-fit parameters for a single Gaussian were 
calculated. 

The XRT and EIS data were aligned by cross correlating the 
Be_thin filter image from XRT with the Fexvi 262.984 A raster 
from EIS. The EIS field of view is shown in the XRT full disk 
image in Figure [1] as solid white lines. The entire EIS field of 
view for the C_poly filter and Fexii 195.119 A raster are also 
shown in Figure [H 

We use magnetic field measurements from the Michelson 
Doppler Imager (MDI) on SoHO. AW of the 96-minute fuh 

disk maPTiptoPTPims t^ikpn within hours of thp H^itpi wprp 
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was selected. This region is shown in Figure [2] with solid lines 
in several example rasters or filter images. The emission in this 
region represents the intensity from the core of the active region 
plus ambient corona between the Sun and the telescope. To re- 
move the ambient corona, we selected four background regions 
at the edge of the field of view, shown in Figure [1] as dashed 
lines. We average the intensity in the region of interest and 
subtract the average background intensity. The background- 
subtracted intensities are given in Tabled EIS intensities are 
given in ergs cm~^ s~^ sr~^ and XRT intensities are given in 
DN s"^ pix"\ 

For this work we have included EIS observations of the Fe XVII 
line at 254.87 A. This line is problematic for several reasons. 
Both the emissivity of this line and the effective area of EIS at 
this wavelength are relatively small, which makes the line dif- 
ficult to observe. Furthermore, the atomic data for this transi- 
tion appears to be inconsistent with the other Fe XVII a nd high 
temperature Ca lines observed with EIS IWarren et al ] (2008 ). 
These differences appear to be about a factor of 2. Since this 
is the highest temperature line observed with EIS outside of 
flares it is useful to accept these limitations so that we can 
have better overlap between the EIS and XRT observations. 

The uncertainties in the EIS intensities given in Table [2] are 
the calculated statistical errors due to photon noise and the 
error in the Gaussian fits and an assumed 22% systematic er- 
ror associated wi th the absolute calibration of the EIS data 
fcang et al. I I2OO6I ). The uncertainties in the XRT intensities 
due to photon noise are more difficult to assess. The method 
used to calculate the XRT errors in Table [2] is described below. 

To generate possible DEM curves that can reproduce the ob- 
served fluxes given in Table[2l we have used xrt_dein_iterative2, 
which was designed originally for use with XRT data only but 
has been modified slightly to allow for inclusion of EIS d ata as 
well . Leading up to the launch of Hinode, Golub et al.l ([2004 ) 
and I Weber et al.l (|2004h tested and validated the method with 
synthetic data. The routine employs a forward- fitting approach 
where a DEM is guessed and folded through each response 
to generate predicted fluxes. This process is iterated to re- 
duce the between the predicted and observed fluxes. The 
DEM function is interpolated using several spline points, which 
are directly manipulated by mpf it, a well-known and much- 
tested IDL routine that performs a Levenberg-Marquardt least- 
squares minimization. This routine uses Monte- Carlo iterations 
to estimate errors on the DEM solution. For each iteration, 
the observed fluxes in each filter were varied randomly within 
the uncertainties and the program was run again with the new 
values. 

The xrt_dein_ inter at ive2 program requires user input for 
the XRT response functions and EIS emissivity functions. The 
XRT filter responses were calculated using the XRT standard 
software (make_xrt_wave_resp and inake_xrt_temp_resp). 
Thpsp nroPTPims ptrrnunt for thp timp Hpnpndpnt rnntam- 
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bake-out was designed to remove the accumulating wavelength- 
dependent contamination layer that was affecting the instru- 
ment sensitivity. As a result, the contamination layer was close 
to maximum thickness. In addition, the XRT team was not yet 
taking regular G-band images, which is one of the tools used 
to estimate the thickness of the contamination layer. As a re- 
sult, the thickness is not as well known. ALmesh observations 
are generally the most affected by contamination, especially at 
low temperatures where the sensitivity is dependent mainly on 
longer wavelength spectral lines. In light of these uncertain- 
ties, and combined with the fact that we have EIS spectral 
lines that effectively cover the lower coronal temperatures, we 
have elected to eliminate ALmesh from the DEM calculation. 

XRT is a broadband instrument allowing photons of many 
different energies to generate electrons on the detector. The 
number of electrons that are generated are proportional to the 
photon energy, so there is no a priori way to deduce the num- 
ber of photons that produced a signal from the number of elec- 
trons deposited onto the detector. The measured count rate in 
a pixel could be from a few high-energy photons or from many 
low-energy photons, and the uncertainty would vary accord- 
ingly. Below we describe how we use a "bootstrap" method to 
calculate the errors associated with photon uncertainties. 

To determine the XRT filter errors due to photon noise, we 
first calculate a DEM using the method described above as- 
suming a 20% error in the XRT filter intensities. We then use 
that DEM to calculate the emerging spectrum from the solar 
plasma, /soZar(A), in photons s~^ cm~^ sr~^ A~^. We convolve 
the spectrum with the effective area for each filter and multiply 
by the resolution, wavelength bin size, and exposure time for 
each filter to find the spectrum in photons at the detector in a 
given filter, i.e., 

Idet{\) = Isolar{\) * EAfilter{\) * AA * texp * C. (l) 

In the above equation, the effective area, EAfuter, is in cm^ 
and is taken from the XRT program make_xrt_wave_resp, the 
AA is the size of the wavelength bin is the exposure 

time, and the factor, C, is the sterradians subtended by one 
XRT pixel. After calculating the photons that interact with the 
detector, we determine the photon noise by taking the square 
root of this intensity, adet — V Idet • The intensity in a given 
pixel in data number (DN) can then be found by 



IdN — ^ IdetW * 



E(\) 1 
eo G 



(2) 



where E is the energy associated with an incoming photon, 
E{X) = /ic/A, Co is the energy per electron on the detector 
(3.65 eV), and G is the gain of the detector (59 electron DN~^). 
We then propagate the photon noise through this step to deter- 
mine the error in the intensity in DN. Note there are additional 
sources of statistical uncertainty in this calculation, but they 



calculated in each spectral line or filter for this DEM are given 
in Table [2l In general, the ratio of the observed to modeled 
intensity is close to 1 with two exceptions. The Fexvil spectral 
line intensity is off by a factor of 2; this is likel y due to poor 
atomic data and is consistent with the results in lWarren et al.l 
(| 2008h . Additionally, the Be_thick intensity calculated from 
the DEM is roughly a factor of 3 higher than observed. 

The dotted lines clustered around the solid thick line are 
DEMs calculated by varying the input intensities within the 
uncertainties and hence provide an estimate of the uncertainty 
of the DEM. Another way to assess the uncertainty in the DEM 
and determine the temperature regime where the DEM is well 
constrained is to calculate how much additional emission could 
be added to a single temperature bin without changing the 
modeled intensities in any spectral line or filter by more than 
the expected errors. This value is shown as a blue line in Fig- 
ure [3l For instance, the DEM calculated for Log T = 7.0 MK 
is approximately 6 x 10^^ cm~^ K~^. We determine we could 
increase the DEM to 2 x 10^° cm"^ K"^ (a factor of 3 x 10^) 
without changing the modeled intensities by more than the er- 
rors. This implies this temperature bin is not well constrained 
and emission could be in this temperature bin, but the spectral 
lines and filters we are using in this analysis are not sensitive 
to it. The ratio of this "maximum DEM" to the calculated 
DEM is less than 3.0 in the temperature range of 6.1 < Log T 
< 6.7; we consider this temperature range well constrained by 
the observations. 

If we define the differential emission measure to be ^ = 
n^ds/dT^ we can write the intensity in a given spectral line or 
filter, /a, in terms of the emissivity function or filter response 
function, ex{T,n), i.e.. 



h 
1 

1 



^ I e,{T,n)^dT 



^ / 
- I 



47r/A 
1 



eioci(T) 



(3) 
(4) 
(5) 



where we have introduced the emission measure loci function. 



6oci(T) 



47r/A 



(6) 



ex{T,n) 

([Jordan et al.lll987l V In the right panel of Figure O the emis- 
sion measure distribution, i.e., ^{T)dT, is shown in black. The 
EM loci curves for the EIS lines considered in this analysis are 
shown in red and the EM loci for the XRT filter intensities are 
shown in green. The blue line shows the maximum emission 
possible in a single bin without changing the modeled intensi- 
ties by more than the errors. Note that the thickest XRT filter 
(Be_thick) does not well constrain the emission measure at high 
temperatures due to the large uncertainties in the intensities. 



Winebarger et al. 



Combining Density Diagnostics and DEM to Constrain Steady Heating 



5 



threshold of 1200 ergs cm ^ s ^ sr ^. The moss regions are 
shown with contours in Figure lU 

Before we calculate the density, we first need to subtract the 
background emission from the moss regions. Because moss 
forms the footpoints of the high temperature loops, we are not 
only looking through the ambient corona, but also through the 
hot loops above the moss. To account for this, we choose to use 
the central region of the core (shown as a rectangle in Figure [2]) 
as the background. These background intensities are the same 
as the apex intensities and given in Tabled 

Using the ratio of the background subtracted intensity and 
the density sensitive emission measure ratio calculated from 
CHIANTI 6.0.1, we calculate a density for each moss pixel. 
The lower left panel shows a density map of the moss regions. 
The density does not appear to depend on spatial location in 
the primary moss regions, though the satellite regions to the 
right typically have lower densities. The lower right panel is 
a histogram of the number of pixels in the moss region with a 
given density. The average density is Log n = 10.15 cm~^ and 
the standard deviation is 0.49. The largest densities measured 
in this region are ^ 5 x 10^° cm~^. 

3.3. Loop Length 

We approximate the loop lengths and geometries in the active 
region core using a potential field extrapolation of the photo- 
spheric field. First, we start with the full Sun, time-averaged 
photospheric magnetic field measurements, shown in the upper 
left panel of Figure [5] We extract a region of full disk mag- 
netic field around the active region, shown with an "X" on the 
figure. The extracted region, shown in the upper right panel 
of Figure \5\ has been generated so that the pixels are approx- 
imately square with the pixel sizes measured in Mm from the 
center of the active region. Hence, we create a magnetic field 
coordinate system that is Cartesian around the middle of the 
active region. We transform the coordinates from the magnetic 
field coordinate system to the image coordinate system using 
the transformation matrixes in Aschwanden et al. (1999). 

After correcting the magnetic field, we numerically solve the 
equations V x B = and V • B = to determine the mag- 
netic field vectors in the volume above the photospheric field 
(e.g., lGarvlll989l V To define the footpoints of the core loops, 
we transform the moss footpoints, originally in the image coor- 
dinate system, to the magnetic field coordinate system. Recall 
that the original moss region was co-aligned with the original 
magnetic field before the region of interest was extracted. Con- 
tours of the moss region in the magnetic field coordinate system 
are shown on the magnetic field image in the top right of Fig- 
ure [5] We then trace field lines from the moss regions. The x,y 
coordinates of the moss are defined by the Fexii 186.880 A in- 
tensity. We assume the height of the moss, and the originating 
z position of the field line, is 2.5 Mm above the photospheric 



of previous successful active region models. 1) We assume 
the geometry of the active region core is well represented by 
field lines from the potential field extrapolation. 2) We as- 
sume the atomic physics of the plasma is well represented by 
the CHIA NTI 6.0.1 package with the ionization balance de- 
scr ibed by Mazzotta et al] (|l998h and abundances described 
bv lFeldman et al.1 (|l99"l i~ 3) We assume the heating along the 
loop is well approximated by uniform, steady heating. 4) We 
assume the cross-sectional area of the loop increases as the 
magnetic field along the loop decrease, i.e., A{s) ^ 1/B(s). 

Below we describe the process of finding the best solution for 
a single field line. We complete this process for each field line, 
then derive a DEM from the portion of all the field lines that 
projects into the region of interest shown in Figure [2] To avoid 
duplicate field lines, we consider field lines traced from only the 
positive polarity magnetic field. 

The example field line is shown as a thick black line in the 
lower left panel of Figure [5] The field line is 36.7 Mm long. 
The (x,z) and (y,z) projections of the field line are shown in 
Figure [6] These figures demonstrate that the field line starts 
and terminates 2.5 Mm above the solar surface and is roughly 
semi-circular in the x-projection. The y-projection shows that 
the field line is slightly inclined. We associate the field line 
with the background subtracted Fexii 186.880 and 195.119 A 
intensities at its originating footpoint, for this field line, 819.5 
and 1213.2 ergs cm~^ s~^ sr~^, respectively. 

The first step in formulating the model is guessing the uniform 
heating rate tha t will ma t ch t he Fexii intensities at the base. 
Recall that Rosn er et al.l (|l97 8) determined scaling laws that 
related the volumetric heating rate to the loop half length in 
cm and base pressure in dyn cm~^, i.e., 

5 7/6 / 7-5/6 



^ = 1 X 10>o 7^ 



(7) 



Using the Fexii line ratio, we calculate a density of 1.0 x 10^° 
cm~^. This density is determined assuming the line intensity is 
generated by plasma at the temperature of peak formation for 
Fexii, Log T = 6.1. We can estimate the base pressure, then, 
by po = 2nokT where k is Boltzmann's constant. For this 
density, we estimate a base pressure of 3.6 dyn cm~^. Using 
this base pressure and the loop half-length in cm, we estimate 
a uniform heating rate of 9.2 x 10~^ ergs cm~^ s~^. 

This value is just an initial estimate of the steady uniform 
heating rate based on assumptions, such as constant loop 
pressure and a simplified radiative loss function. Using this 
estimate, we solve the one-dimensional hydrodynamic equa- 
tions of co ntinuity, momentum, and energy with a steady- 
state solver (jSchriiver &: v an Ballegooiien 2005) . The radiative 
loss function we use to solve the equations was described by 
brooks WarrenI (|2006h . The density and temperature along 
the loop for this heating rate is shown in the bottom two panels 
in Figure [6] as solid lines . 
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of 4910 ergs cm ^ s ^ sr ^ and a Fexii 195.119 intensity of 
6950 ergs cm~^ s~^ sr~^. 

The intensity found in Equation [8] would be the intensity if the 
plasma filled the area of the pixel. To determine the actual per- 
centage of the pixel that contains plasma, or the filling factor, 
we take the ratio of the observed to modeled Fexii 195.119 A 
intensities. For this simulation, the filling factor is 0.178. Then 
we multiply the simulated Fexii 186.880 A intensity by this fill- 
ing factor. The Fexil 186 A intensity becomes 873 ergs cm~^ 
s~^ sr~^ which is slightly larger than the observed 819.5 ergs 
cm~^ s""*^ sr""*^. 

Because the simulated intensity in Fexii 186.880 A is too 
large, we reduce the heating rate by a factor of 2; if it had 
been too small, we would have increased the heating rate by a 
factor of 2. Using the new heating rate, we solve again the one- 
dimensional hydrodynamic equations. The temperature and 
density solution for a heating rate of 4.6 x 10~^ ergs cm""^ s~^ 
is shown in the lower panels of Figure [6] with dashed lines. This 
solution produced a Fexil 186.880 A intensity of 755 ergs cm~^ 
sr-1. 

We now have two different heating rates and two differ- 
ent Fexii 186.880 A intensities that bracket the desired Fexii 
186.880 A intensity. We interpolate between them and predict 
a heating rate of 6.8 x 10~^ ergs cm""^ s~^ will achieve the de- 
sired Fexii 186.880 A intensity. The temperature and density 
solution for this heating rate is shown in Figure [6] with dash-dot 
lines. The Fexil 186.880 and 195.119 A intensities determined 
from Equation [8] were 1791 and 2702 ergs cm~^ s~^ sr~^, re- 
spectively. Using the Fexii 195.119 A ratio, we determine a 
filling factor of 0.457. The simulated Fexii 186.880 A intensity 
is then 819 ergs cm~^ s~^ sr~^ which matched the observed 
intensity. 

We complete this process for all field lines in our study. For 
10% of the field lines, the necessary filling factor was > 1; we 
eliminated these field lines from the study. After the best solu- 
tion was found for all the field lines, a DEM of was generated 
from the portion of the density and temperature that projected 
into the region of interest (the white rectangle shown in Fig- 
ure [2]) Figure [71 shows the comparison of the DEM calculated 
from the data (black) and the DEM predicted by the steady 
heating model (green). The model DEM well approximates the 
observed DEM around its peak, but the model DEM does not 
well describe the observations at lower or higher temperatures. 
We discuss this discrepancy in the next section. 

Through this process, we solved for the density and temper- 
ature along the loop, as well as the required filling factor to 
bring the Fexii intensities into agreement with observations. 
Figure [8] shows the resulting relationship between the filling 
factor and temperature (left panel) and loop length and tem- 
perature (right panel). 



equations for steady heating, we adjust the magnitude of the 
heating rate iteratively until the simulated intensities in the 
Fexii 186.880 and 195.119 spectral lines match the observed 
intensities to within 1%. We then construct a DEM from the 
solutions and compare with the observed DEM. 

In the temperature range 6.3 < Log T < 6.7, the model 
DEM is in general agreement with the observations. The mag- 
nitude of the model DEM is approximately the same as the 
magnitude of the observed DEM. The median temperature of 
the model DEM (Log T = 6.5 MK) is the same as the tem- 
perature of the peak of the observed DEM. This good agree- 
ment was arrived at using relatively few assumptions: poten- 
tial field geometry, CHIANTI atomic physics, uniform, steady 
heating, and loop area inversely proportional to the magnetic 
field. These assumptions were motivated by successful param- 
eter space searches in previous studies. 

Unlike previous studies, however, no additional "fudge fac- 
tors" were used, such as a filling factor, to force agreement 
between the observations and the model. Because this study 
relies on density sensitive line intensities, the areal filling fac- 
tor is determined through the model calculation. After making 
these four key assumptions, there are no additional assump- 
tions that could be made that could, for instance, artificially 
raise or lower the temperatures in the loops. In fact, the max- 
imum and minimum temperature of model DEM are defined 
by the maximum and minimum density and loop length shown 
in Figures [H and [5] We determined the maximum electron 
density in the moss was 5 x 10^° cm"*^ and the maximum 
loop length measured from potential field extrapolations was 
^ 100 Mm. The moss densities were measured with density 
sensitive Fe Xll lines formed at 1.5 MK. For the resulting pres- 
sure (^ 20 dyn cm^) and half-length, we estimate the maximum 
apex temperature of loops for high-frequency heating is 6.5 MK 
from the RTV scaling law which is in agreement with the max- 
imum temperature of the model DEM. Observations of signif- 
icant emission at higher temperatures would have indicated a 
disagreement with the steady heating model. 

At temperatures much lower or higher than the peak tem- 
perature (Log T < 6.1 or Log T > 6.7), the observed DEM is 
not well constrained and comparisons with the model are not 
meaningful. At mid-range, "warm" temperatures (6.1 < Log 
T < 6.3), however, the observed DEM is well constrained and 
significantly larger than the DEM from the model. There are 
(at least) two possible explanations for this discrepancy. 1) 
The warm emission is from the overlying arcade of warm loops 
that was not removed in the original background subtraction. 
Figure [T] supports this view. The Fe xii raster shows the pres- 
ence of many overlying loops that are not confined to the core 
loops bright in the XRT C_poly image. 2) The warm emission 
is truly from the core loops, casting serious doubt on the steady 
heating model. Another heating model , such as the infrequent 
heating or "nanoflare' model (Carsill Klimchuk 1997) must 
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ture emission at the apex of the core loops, and the appearance 
of warm loops in the series of TRACE images. Regardless, ad- 
ditional modeling efforts are underway to t est whether the low- 
frequency nanoflare model (as described by Cargi ll KlimchukI 
fT?97, 2004; Khmchuk et al. 2008; Khmchuk 2009^ could be op- 
erating in active region cores. Initial studies characterizing 
the DEM p redicted b y long nanoflare storms are underway 
(iKhmchuk et al. 20 081: ISusino et al.ll2Q10l: iTripathi et al.ll201li : 
IWarren et al.ll2011l : iMuTu-Moore et al.ll201ll ~ 

In this paper, we chose to use the cooler Fexii 186.880/195.119 A 
line ratio in the moss region and the hotter Fe xiil 202.044/203.826 . 
line ratio to determine the density in the core region. It 
has been well documented, however, that densities calculated 
from the Fexii line ratio could be as much as a factor of 2 
more than the densities calculated from other EIS line ratios 
(| Warren et al.|[2010l l. For steady heating, the apex temper- 
ature of a loop is proportional to the cube root of the base 
pressure, i.e., Tapex ^ (po)^^^, so if the density were a factor of 
2 lower, the temperatures in the loops would be a factor of 1.25 
lower which would shift the model DEM to lower temperatures. 
Currently, the average median temperature of the simulated 
loops is 3.1 MK. If the true densities were a factor of 2 lower, 
it would shift this average temperature to 2.5 MK. 

The densities measured in this analysis are in go od agreement 
with other measurements of moss densities (Tripathi et al. 
I2OO8I , [ioiol V The emission measure curve found in this analysis, 
however, is significantly steeper in the temperature range of 6.0 
< Log T < 6.5 than have been found in previous emission mea- 
sure calculations. We find the emission measure curve in this 
temperature range can be approximated as a power-law, i.e., 
EM - T^ with an i ndex, 6, of 3.2, while p revious results have 
found 1 < 6 < 3 fe.g.. lDere MasonllT993l : lBrosius et al.lll996l ). 
These previous results averaged the intensities over large fields 
of view which included both high temperature loops, moss, 
and extended EUV loops. Two recent analyses of the emission 
measure distribution of int er-moss regions determined a range 
of in dices 2.1 < b < 3.4 (|Warren et al.1 12OIII : iTripathi etall 
l201lh . It is clear a systematic study of the emission measure 
distributions of inter-moss regions is required to fully charac- 
terize the temperature structure of the active region core. 

In this paper, we combine XRT and EIS data s ets to cal- 
culate a differential emission measure, similar to Sch melz et al.l 
(|2010l V Figure [9] demonstrates the power of this combined data 
set compared to using the data sets individually. The DEM 
shown in blue is calculated from the XRT intensities alone, the 
DEM shown in red is calculated from the EIS line intensities 
alone, and the DEM shown in black is the combined data set. 
Comparing the combined DEM with the individual instrument 
DEMs, we see the combined DEM agrees well with the EIS 
DEM in the range 6.0 < Log T < 6.5, however the EIS DEM 
greatly overestimates the true DEM at high temperatures. The 
DEM calculated from the XRT data is significantly different 



perature overlap of the XRT filters. The data set presented 
in this analysis, however, contained only the Fexvii 254.870 A 
line that overlaps the XRT temperature range. Furthermore, 
additional studies of combi ned dat a sets that did ha ve signif- 
icant temperature overlap (|Warren et al. 2 oTll : lO'Dwv er et al] 
l201lh did not find the need to include a cross-calibration factor. 
Because the XRT contamination varies in time, it is diflficult to 
assess the relative calibration of the two instruments. 

Finally, our DEM curve for AR 10955, which used XRT data 
from 16:00 UT (see Table 1), shows no significant hot plasma 
with Log T > 7.0; the blue line in Figure 3 shows our limited 
sensitivity to such plasma, given the observed count rates and 
errors in the count rates. This active regio n, however, has been 
studied previously. I Schmelz et al ] (|2009bl ) used a similar XRT 
data set from 18:00 UT (approximately two hours after the 
XRT data considered in this paper) to determine the DEM for 
a port ion of t he active region to the north-west of the core. 
Schmel z et aD (^09a) used the 18:00 UT XRT data as wefl 
as RHESSI upper limits to modify and further constrain the 
DEM. Both papers found that considerable emission in the high 
temperature range was required to account for the small but 
significant signal detected in the XRT Be_thick filter. 

Although both the 16:18 UT and 18:00 UT Be_thick images 
had the same resolution, field-of-view, and exposure time, the 
earlier image analyzed in this paper does not show a significant 
signal, where the later one analyzed by Schme lz et a l. (2009b, a) 
does. The GOES signal for 2007 May 13 was low ah day, rarely 
getting above level AO. There were two small A4 fiares from 
AR 10955 at 11:20 and 11:35 UT, before both sets of XRT 
observations, but then the region settled down to its sub- AO 
level again. It maintained the quiescent stage through our 16:00 
UT observations, only to rise again to A1-A2 from 17:00 UT to 
just past 18:00 UT when there was another small A7 fiare. One 
possible explanation for the different DEM results is that AR 
10955 appears to have been in a more quiescent state during 
our 16:00 UT observations, and in a somewhat heightened state 
of activity (right before the small A7 fiare) during their 18:00 
UT observations. It would be interesting to repeat the XRT- 
RHESSI DEM analysis for a stronger, more powerful active 
region with a significant signal in the XRT Be_thick filter to 
see if the high-temperature plasma is detectable. 

Hinode is a Japanese mission developed and launched by 
ISAS/JAXA, with NAOJ as domestic partner and NASA and 
STFC (UK) as international partners. It is operated by these 
agencies in co-operation with ESA and the NSC (Norway). 

ARW was supported by an NSF Career grant. Solar physics 
research at the University of Memphis is supported by a Hin- 
ode subcontract from NASA/SAO. ARW thanks Mark Weber 
for many enlightening conversations on Differential Emission 
Measures. 
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Log T (K) Log T (K) 

Fig. 3. — (Left Panel:) The thick sohd hnes show the best DEM for the input XRT and EIS intensities and the dotted hnes 
surrounding the sohd hne gives an estimate of the error in the DEM. The blue hne represents the maximum emission that can be 
added to a single temperature bin without changing the modeled intensity by more than the observed errors. {Right Panel:) The 
integral form of the DEM and corresponding EM loci curves. The EIS EM loci curves are shown in pink and the XRT EM loci 
curves are shown in green. 



Table 2 
EIS AND XRT Intensities 



Line/Filter 
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Idem 


I/Idem 


Fex 177.239 


254 


194 


239 


1.1 


Fexi 180.401 


482 


218 


637 


0.76 


Fexii 186.880 


587 


149 


434 


1.4 


Fexii 195.119 


768 


259 


941 


0.81 


Fexiii 203.826 


1725 


414 


1153 


1.5 


Fexiii 202.044 


601 


213 


421 


1.4 


Fexiv 264.787 


1016 


241 


1189 


0.85 


Fexiv 274.203 


870 


215 


864 


1.0 


Fexv 284.160 


9450 


2150 


11570 


0.82 


Fexvi 262.984 


1140 


257 


852 


1.3 


Fexvii 254.870 


33 


16 


58 


0.57 


ALmesh'' 


1940 


390 


1406 


1.4 


C-poly 


1226 


247 


1075 


1.1 
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Fig. 7. — DEM calculated from XRT and EIS data (black). DEM calculated from a steady heating model (green). 
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Fig. 9. — DEM calculated from XRT intensities only (blue), from EIS intensities only (red) and from the combined XRT and EIS 
data set (black). 



